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Kaon decays have played a key role in the construction of the Standard Model and continue to 
be an important testing ground of the electroweak flavour theory. They can provide new signals 
of CP-violation phenomena and, perhaps, a window into physics beyond the Standard Model. 



The interplay of long-distance QCD effects in strangeness-changing transitions can be analyzed, 
combining the short-distance Operator Product Expansion with Chiral Perturbation Theory tech- 
niques. A brief overview is presented, focusing on a few selected decay modes. A more detailed 
r ^ and comprehensive review can be found in Ref. [1], 
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1. Theoretical Framework 

Kaons have been at the center of many fundamental developments in particle physics, playing 
a key role in the construction of what we now call the Standard Model (SM): the introduction of 
internal flavour quantum numbers (strangeness) [2, 3], parity violation [4, 5], meson-antimeson 
mixing [6, 7], quark mixing [8, 9], the discovery of CP violation [10], and the suppression of 
flavour-changing neutral currents and the GIM mechanism [11]. High-precision experiments on 
rare kaon decays provide sensitivity to short-distance scales (c, t, W ± , Z) and offer the exciting 
possibility of unravelling new physics beyond the SM. Searching for forbidden flavour-changing 
processes beyond the 1(T 10 level [Br(# L -> e ± fi^) < 4.7 x 1(T 12 [12], Br(K L -> e ± e ± n T n T ) < 
4.12 x 1(T 11 [13], Br(K+ -> K + jX + e-) < 1.3 x KT 11 [14], Kr{K + -> 7i+p-e+) < 5.2 x 1(T 10 [15] 
(90% C.L.)], one is actually exploring energy scales above the 10 TeV region. The study of allowed 
(but highly suppressed) decay modes provides, at the same time, very interesting tests of the SM 
itself. Electromagnetic-induced non-leptonic weak transitions and higher-order weak processes 
are a useful tool to improve our understanding of the interplay among electromagnetic, weak and 
strong interactions. In addition, new signals of CP violation, which would help to elucidate the 
source of CP- violating phenomena, can be looked for. 

The theoretical analysis of kaon decays is highly non-trivial. While the underlying flavour- 
changing transitions among the constituent quarks are associated with the electroweak scale, the 
corresponding hadronic amplitudes are governed by the long-distance behaviour of the strong in- 
teractions, i.e., the confinement regime of QCD. The short-distance approach to weak transitions 
makes use of the asymptotic freedom property of QCD to successively integrate out the fields 
with heavy masses down to scales /I < m c . Using the operator product expansion (OPE) and 
renormalization-group techniques, one gets an effective AS = 1 Hamiltonian [16] 

Kf =x = ^v^viEGooa+hx., Clip) = soo-^oo^s ( i.d 

V2 t V ud V us 

which is a sum of local four-fermion operators <2r, constructed with the light degrees of free- 
dom (u,d,s;e,/J.,Vi), modulated by Wilson coefficients C,-(ju) which are functions of the heavy 
(Z,W,t,b,c,x) masses. The CP-violating decay amplitudes are proportional to the components 
yi(H). The overall renormalization scale /I separates the short- (M > /i) and long-distance (m < /i) 
contributions, which are contained in Q(/i) and <2;, respectively. The Wilson coefficients are 
fully known at the next-to-leading order (NLO) [17, 18]; this includes all corrections of &(a%t n ) 
and tff(a" +1 t"), where t = \og{M\/M2) refers to the logarithm of any ratio of heavy mass scales 
(M\ t 2 > At)- In order to calculate the kaon decay amplitudes, we also need to know the nonpertur- 
bative matrix elements of the operators <2, between the initial and final states. 

The low-energy strong interactions are better understood with symmetry considerations, be- 
cause the lightest pseudoscalar mesons correspond to the octet of Goldstone bosons associated with 
the dynamical chiral symmetry breaking of QCD: SU (3)l 8) SU (3)r — > SU (3)y. Their dynamical 
properties can then be worked out systematically through an effective Lagrangian. The quark and 
gluon fields of QCD are replaced by a unitary matrix U((f>) = exp(i\/2<i> / F ), parameterizing the 
Goldstone excitacions over the vacuum quark condensate {q J L q' R } = u,d,s). The Chiral Per- 
turbation Theory [19, 20, 21, 22] (#PT) formulation of the Standard Model is given by the most 
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Figure 2: Short-distance Feynman diagrams. 





Figure 1: Evolution from M\y to the kaon mass scale. Figure 3: Long-distance Feynman diagrams. 



general effective Lagrangian, involving the matrix U(ty), which is consistent with chiral symmetry. 
The Lagrangian can be organized in terms of increasing powers of momenta (derivatives) and quark 
masses over the chiral symmetry-breaking scale (A^ ~ 1 GeV). All short-distance information is 
encoded in the low-energy couplings (LECs) of the #PT operators. At lowest order (LO), &(p 2 ), 
the strong interactions are fully parameterized in terms of only two LECs: the pion decay constant 
F ~ F n = 92 A MeV and Bo — — < uu > /F 2 , which accounts for the explicit chiral symmetry 
breaking through the quark masses. Ten additional LECs L,- are needed at &(p 4 ). 

Non-leptonic weak interactions with AS = 1 are incorporated as a perturbation to the strong 
#PT Lagrangian. At LO the most general effective Lagrangian, with the same SU (3)l <8> SU (3)r 
transformation properties as the short-distance Hamiltonian (1.1), contains three terms [1]: 

^f =l = ~^V ud V: s |g 8 (XL,L») + g 21 + h, 2l L^ + e 2 g &gew F 6 (W'QU) +h.cj , 

(1.2) 

where = iF 2 U^D^U represents the octet of V — A currents, A = — iXi)/2 projects onto 
the s — » d transition, Q = ^diag(2, — 1, — 1) is the quark charge matrix and (• • •) denotes the 3- 
dimensional flavour trace. The LECs and g 2 j measure the strength of the two parts of J^gr 
transforming as (8/., Ir) and (27^, 1#), respectively, under chiral rotations, while g ev/ accounts for 
the electromagnetic penguin operators. 

The #PT framework determines the most general form of the K decay amplitudes, compati- 
ble with chiral symmetry, in terms of the LECs multiplying the relevant chiral operators. These 
LECs, which encode the short-distance dynamics, can be determined phenomenologically and/or 
calculated in the limit of a large number of QCD colours (matching). Chiral loops generate non- 
polynomial contributions, with logarithms and threshold factors as required by unitarity. Fig. 1 
shows schematically the procedure used to evolve down from Mw to ra# . While the OPE resums 
the short-distance logarithmic corrections log (M/jj.), the #PT loops take care of the large infrared 
logarithms log(/x jm % ) associated with unitarity corrections (final-state interactions). 
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\gjijge\ Source 

1.0021 ±0.0016 K-^n/e 
0.9978 ±0.001 8 K->n/e 
1.0010 ±0.0025 K-+np/e 
1.0018±0.0014 T->n/e 



Table 1: Determinations of the ratio of the jx 
and e~ couplings to the W. 
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Figure 4: Constraints on \V uc t\ and |V„ S | from Ka de- 
cays (horizontal), + —> + nuclear decays (vertical) 
and K^ajit^a (oblique band) [33]. 



2. Leptonic and Semileptonic Decays. 

Leptonic and semileptonic kaon decays are well understood theoretically, including electro- 
magnetic corrections. Strong interactions only appear through the hadronic matrix elements of the 
left-handed weak current, which can be precisely studied within #PT and with lattice simulations. 

The ratios R^) = T[P e~V e (y)]/T[P~ -> p^v^y)} (P = n,K) have been calculated 
[23, 24] and measured [25, 26, 27, 28, 29] with high accuracy, allowing for a test of charged- 
current lepton universality at the 0.2% level. As shown in Table 1, similar precisions have been 
achieved in K — > nlvg [30] and x — > v r £ve [31] decays. 

The ratio of radiative inclusive decay rates F[K — > /i _ v iU (7)]/r[7r _ — > pL^v^y)] provides 
information on the quark mixing matrix [30, 32]. With a careful treatment ot electromagnetic and 
isospin-violating corrections, one extracts IV^/V^I \Fk/F% \ = 0.2763 ±0.0005 [33]. Taking for the 
ratio of meson decay constants the lattice average Fk/F k = 1.193 ±0.006 [34], this gives 

\V US /V ud \ = 0.2316 ±0.0012. (2.1) 

The most recent Kg3 experiments have resulted in improved precision and significant shifts in 
the branching fractions [36]. Including electromagnetic and isospin-breaking corrections [37, 38], 
one obtains \V m f + (0)\ = 0.2163 ± 0.0005 [30], with / + (0) = 1 ± 0[(m s - m u ) 2 } the K° -> K-£ + v e 
vector form factor. The exact value of / + (0) has been thoroughly investigated since the first precise 
estimate by Leutwyler and Roos, /+(0) = 0.961 ±0.008 [39]. While analytical calculations based 
on ^PT obtain higher values [40, 41], owing to the large (~ 0.01) 2-loop corrections [42], lattice 
results [34] tend to agree with the Leutwyler-Roos estimate. Taking as reference value the most 
recent and precise lattice result [43], /+(0) = 0.960 ±0.006, one obtains [1] 

\Vus\ = 0.2255 ± 0.0005 exp ± 0.0012 th . (2.2) 

Together with \V ud \ = 0.97425 ± 0.00022 [35] and the negligible \V ub \ contribution [36], the deter- 
minations (2.1) and (2.2) imply an stringent test of the unitarity of the quark mixing matrix [33]: 

A C KM = \V ud \ 2 + \V lls \ 2 + \V ub \ 2 -l = 0.0001 ±0.0006. (2.3) 
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3. Nonleptonic Decays and Direct CP Violation: e'/e 

The measured A (K — > nn)i decay amplitudes show a strong enhancement of the octet AI = A 
transition amplitude into a 2% final state with isospin 1 = 0: IA0/A2I ~ 22. In the ^PT framework 
this manifests as a huge difference between the LECs in Eq. (1.2). A LO fit to the data gives 
\g$\ ~ 5.0 and |g27| — 0.285. Part of the enhancement originates in the strong rescattering of the 
final pions, which at one loop increases Ao by roughly 35%; taking the ^PT 1-loop contributions 
into account, one finds a sizeably smaller octet coupling (the central values change slightly to 3.61 
and 0.297, respectively, if isospin violation is included) [1, 44, 45]: 

\g&\ = 3.62±0.28, \g 27 \ = 0.286 ±0.28. (3.1) 

In the absence of QCD corrections, the SM (W exchange) implies g& = gn = |, very far from 
the phenomenologically required values. The computed short-distance QCD corrections show the 
needed qualitative trend to understand the data, but a proper calculation of the hadronic matrix 
elements of the relevant Qj operators is still lacking. The matching of the effective descriptions 
^f/ =1 (short-distance) and Jzf 2 A5=1 (#PT) can be done in the large-Afc limit with the result g% = 
1.13 ±0.18 and gjj = 0.46 ±0.01 [1], which shows the relevance of the missing NLO corrections 
in 1 /Nc- Lattice simulations have recently achieved a quite successful description of A2 (gn) [46], 
but a real quantitative understanding of Ao (g%) remains still problematic [47]. 

The situation is much better fot the ratio [48, 49, 50, 51] [r\ ab = A (Kl — > % a Ti b )/A(K s -> % a n h )\ 



Re (e'/e) = l -[ 



16.8 ±2.0) x 10" 4 , (3.2) 



which demonstrates the existence of direct CP violation in the K — > 2% decay amplitudes. When 
CP violation is turned on, the amplitudes A/ acquire imaginary parts. To first order in CP violation, 

ImAo ImA2 



V2 ReA c 



ReA ReA 2 



(3.3) 



where the strong phases %i can be identified with the S-wave Tin scattering phase shifts at = tuk , 
up to isospin-breaking effects [44, 45]. The phase (j> E i =%2 — %o + ft/2 = (42.5 ±0.9)° is very close 
to the so-called superweak phase, <f> e w tan -1 [2(rriK L — rriK s ) / (r^ s — Tk l )] = (43.51 ± 0.05)°, im- 
plying that cos ((j> e i — (j> e ) w 1. The CP-conserving amplitudes ReA/ can be set to their experimen- 
tally determined values, avoiding in this way the large uncertainties associated with the hadronic 
matrix elements of the four-quark operators in . Thus, one only needs a first-principle 

calculation of the CP-odd amplitudes ImAo and ImA2; the first one is completely dominated by 
the strong penguin operator while the leading contribution to the second one comes from the 
electromagnetic penguin Q%. Fortunately, those are precisely the only operators that are well ap- 
proximated through a large-A^ estimate of LECs, because their anomalous dimensions are leading 
in 1 /Nc- Owing to the large ratio ReAo /ReA2, isospin violation plays also an important role in e' / £ 
[45]. The one-loop #PT enhancement of the isoscalar amplitude [52, 53] destroys an accidental 
LO cancellation of the two terms in (3.3) [54, 55, 56], bringing the SM prediction of e'/e in good 
agreement with the experimental measurement in Eq. (3.2) [52, 53, 57]: 

Re (e'/e) = (\9±2^ 9 6m ±6 1/JVc ) x lO" 4 . (3.4) 



5 



Kaon Physics 



Antonio Pich 




Figure 5: Lowest-order contributions to K$ — > yy. 

4. Rare and Radiative Decays 

Kaon decays mediated by flavour-changing neutral currents are suppressed in the SM and their 
main interest, other than their own understanding, relies on the possible observation of new physics 
effects. Most of these processes are dominated by long-distance contributions; however, there are 
also processes governed by short-distance amplitudes, such as K — > nvv. 

4.1 K° -»• yy and K° -> t+f 

The symmetry constraints do not allow any tree-level K®yy coupling at &{p A ) (X{ } 2 are the CP- 
even and CP-odd states). The decay Ks — > yy proceeds then through a one-loop amplitude, with 
intermediate 7t + 7l~, which is necessarily finite because there are no counterterms to renormalize 
divergences. The resulting 0{p A ) prediction [58, 59], Br(K$ — > yy) = 2.0 x 10 , is slightly lower 
than the experimental measurement Br(A"s — > yy) = (2.63 ±0.17) x 10~ 6 [36]. Full agreement is 
obtained at &(p 6 ), once rescattering corrections (Ks — 5- 7T7T — > n + n~ — > yy) are included [60]. 

The 2-loop amplitude K$ — > y*y* — > £ + £~ is also finite [61] because chiral symmetry forbids 
any CP-invariant local contribution at this order. The predicted rates, Br(Ks — > e + e~ ) = 2. 1 x 10 14 
and Br(^5 — > jU + /i~) = 5.1 x 10~ 12 [61], are well below the experimental upper bounds Br(^s — > 
e + e~) < 9 x 10~ 9 and Br^ -)• /i+aO < 3.2 x 10~ 7 (90% C.L.) [36]. This calculation allows 
us to compute the longitudinal polarization Pi of either muon in the decay Kl — > jU + jU~, a CP- 
violating observable which in the SM is dominated by indirect CP violation from K°-K° mixing. 
One finds \P L \ = (2.6 ±0.4) x 10~ 3 [61]. 

4.2 K -)■ Jiyy 

Again, the symmetry constraints do not allow any tree-level contribution to K2 — > K yy from 
G{p A ) terms in the Lagrangian. The decay amplitude is therefore determined by a finite loop cal- 
culation [62, 63, 64]. Due to the large absorptive % + %~ contribution, the spectrum in the invariant 
mass of the two photons is predicted to have a very characteristic behaviour (dotted line in Fig. 7), 
peaked at high values of m rr . The agreement with the measured distribution [65] is remarkably 
good. However, the @{p A ) prediction for the rate, BR(^ L — > n°yy) = 6.8 x 10~ 7 [62], is signifi- 
cantly smaller than the present PDG average [36] 

BR(K L ^7t°yy) = (1.27 ±0.03) x 10" 6 , (4.1) 
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indicating that higher-order corrections are sizeable. Unitarity corrections from Kl K + K~K 
[66, 67] and local vector-exchange contributions [66, 68] restore the agreement at <^(/? 6 ). 

A quite similar spectrum is predicted [69] for the charged mode — > n ± yy, but in this case 
there is a free LEC already at 0(p 4 ). &(p^) corrections have been also investigated [70]. Both the 
spectrum and the rate can be correctly reproduced [71, 72]. 

4.3 K L -> K°e + e- 

This decay is an interesting process in looking for new CP- violating signatures, because — > 
7i°y* violates CP [69, 73]. The CP-conserving amplitude proceeds through a 2y intermediate state 
and is suppressed by an additional power of a. Using the Kl — > 7i°yy data, the CP-conserving 
rate is found to be below 10~ 12 [1]. The Kl — > n°e + e~ transition is then dominated by the 0(a) 
CP- violating contributions [69], both from K°-K° mixing and direct CP violation. The estimated 
rate Br(^ L -> n°e + e-) = (3.1 ±0.9) x 10~ n [1, 74, 75] is only a factor 10 smaller than the present 
(90% C.L.) upper bound of 2.8 x 10~ 10 [76] and should be reachable in the near future. 

4.4 K — > 7TVV 

Long-distance effects play a negligible role in K + — > 7t + vv and Kl — > 7i°vv. These processes 
are dominated by short-distance loops (Z penguin, W box) involving the heavy top quark; the 
charged mode receives also sizeable contributions from internal charm-quark exchanges. The decay 
amplitudes are proportional to the hadronic matrix element of the AS = 1 vector current, which 
(assuming isospin symmetry) can be obtained from Ka decays: 

T(K -> ttvv) ~ £ F(y id V*; Xi ) (v L y^v L ) (K\s L fd L \K) , Xi = n%/M& . (4.2) 

i=c,t 

The small long-distance and isospin- violating corrections can be estimated within #PT. The Kl — > 
7i°vv transition violates CP and is completely dominated by direct CP violation, the contribution 
from K°-K l) mixing being only of the order of 1%. Taking the CKM inputs from global fits, one 
predicts Br(K L -»• 7T°vv) = (2.4 ±0.4) x 10" 11 and Br(#+ 7T+vv) = (0.78 ±0.08) x 10" 10 
[77, 78]. The uncertainties are largely parametrical, due to CKM input, m c , m t and a s (Mz). 

The K + mode has been observed [79] , while only an upper bound exists on the Kl one [80] : 

Br(K + ^tt+vv) = (1.73+J;^) x lO" 10 , Br(£ L -> 7T°vv) < 2.6 x 10~ 8 (90%C.L.). (4.3) 
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New experiments, aiming to reach £?(100) events (assuming SM rates), are under development 
at CERN (NA62) and J-PARC (KOTO) for charged and neutral modes, respectively. Increased 
sensitivities could be obtained through the recent ORKA proposal for a K + — > 7T + VV experiment 
at Fermilab and the higher kaon fluxes available at Project-X [81]. 
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